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Using a published atlas of reflectance curves and physicochemical 
properties of soils, a statistical analysis was carried out. Reflectance 
bands which correspond to five of the wavebands used by NASA's Thematic 
Mapper were examined for relationships to specific soil properties. The 
properties considered in this study include: Sand Content, Silt Content, 

Clay Content, Organic Matter Content, Cation Exchange Capacity, Iron 
Oxide Content and Moisture Content. Regression of these seven properties 
on the mean values of five TM bands produced results that indicate that 
the predictability of the properties can be increased by stratifying the 
data. The data was stratified by parent material, taxonomic order, 
temperature zone, moisture zone and climate (combined temperature and 
moisture). The best results were obtained when the sample was examined 
by climatic classes. The middle Infra-red bands, 5 and 7, as well as the 
visible bands, 2 and 3, are significant in the model. The near Infra-red 
band, band 4, is almost as useful and should be included in any studies. 
General linear modeling procedures examined relationships of the seven 
properties with certain wavebands in the stratified samples. These 
results reinforced the hypothesis that the TM bands can be utilized to 
predict certain soil properties. Some relationships between curve shape 
and soil properties were also investigated and produced positive results. 
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INTRODUCTION 


Spectral reflectance of soil is influenced by the physical and chemical 
properties of that soil. Some of these physiochemical properties are more 
influential than others. These properties can be measured in the laboratory 
using standard tests. Recent works [Stoner and Baumgardner, 1980 ], [Crouse, 
et al., 1983] i [Kiesewetter , 1982], have suggested relationships of several 
properties to the spectral signature of the respective soils, particularly 
in certain wavebands. The objective of these studies has been to utilize 
the Landsat Satellite data in ’identifying and classifying soils of the 
earth. Stoner ( 1979 ) measured the spectral reflectance of 485 soils of the 
United States along with physical, chemcial and engineering properties of 
those soils. The result of that investigation has been published as the 
"Atlas of Soil Reflectance Properties," [Stoner, et a*.., 1979 ]. This data 
was used by Kiesewetter ( 1982 ) and Crouse, et al. ( 1983 ) » to statistically 
establish relationships between the soil properties and reflectance at 
certain wavelengths. Crouse has selected six wavebands which correspond to 
those of the Thematic Mapper (TM) aboard NASA's Landsat 4 Satellite. That 
study suggested that four, and possibly 5 , of the six TM wavebands may be 
important in the identification of certain properties. Figure shows a 
typical soil presentation from Stoner's work. The properties which showed 
the best results in regression studies were, in order of their relevance, as 
follow: Organic Carbon Content, Water Content, Clay Content, Cation Exchange 

Capacity (C.E.C.), Calcium, Magnesium, Fine Sand, Case Saturation, Sand 
Content, Extractable Acidity, and Iron Oxide Content. In addition, ten 
other properties were included in the original regression analysis but were 
later dropped from consideration because they showed little potential of 
being predicted from TM data. These properties include: Potassium, Medium 


Sand Content, .pa Coarse Sand Content, Silt Content, Aluminum,, Coarse Silt, 
Fine silt, Very Coarse Silt, Sodium; ana Very Fine Sand. 

This current work looks at several properties and five of the six TM 
^ wavebands which are in the range of spectral reflectance measured by Stoner 

5 (1979) and attempts to stratify the soils on the basis of Climate, Taxonomy 

and Parent Material. The objective of this investigation is to determine 
which soil properties show an ability to be predicted from T.M., spectral 
measurements, knowing some information about me temperature and moisture 
zor.es (climate) and the underlying pare .c ma«.'rial. This information could 
conceivably be available for an unkno m soil from sources such as a geographic 


atlas or geologic maps and reports. 



THEORY 


The physical and chemical properties of the soil, in combination, 
produce the spectral reflectance. Some of these properties exert more 
influence than others. Previous work, (Stoner and Baumgardner, 1980) sug- 
gests that moisture content, organic matter, iron oxide content, texture, 
cation exchange capacity are more influential than other properties. 
Kiesewetter (1982) demonstrated that strong negative correlations exist when 
organic matter, iron oxide, and moisture content are compared with reflectance 
values. As the percentage of each of these increases the reflectance 
decreases. The strongest negative correlation seems to be with moisture 
content. Crouse et. al., (1983) have found, through regression analysis, 
that eleven properties seem to be dominant in determining reflectance. 

Since reflectance is dependent upon the physical and chemical makeup of 
the soil, it is logical to try to ascertain how these properties are acquired. 
Parent material should be a place to begin since the soil will inherit 
certain characteristics from this source material. We know, however, that 
soils from similar parent materials can vary greatly in physical and chemical 
composition so other factors must be in play at the same time. Climate 
(temperature and moisture) is a contributing factor since temperature and 
moisture influence biological and chemical processes at work in a soil. 

Parent material is significant in determining the soil texture and chemical 
composition. Parent material may be residual (weathered from underlying 
rock or sediments) or it may be transported by means of wind, water or ice 
and deposited on the bedrock. Other factors at work in the formation of 
soils are topography, living organisms and time. Time is an important 
factor since it changes the characteristics of weathered material from those 
which are inherited from the parent material to those which are acquired in 
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the aging process. The presence of organisms, both plant and animal, is a 
strong contributor to the overall reflectance of the soil. 

Another consideration was to examine the taxonomic order of the soil. 

This is done, primarily, to investigate the possibility that soils classi- 
fied in a particular order will have similar reflective characteristics. 
Taxonomic order, as defined by the DSDA in 1972, is a classification system 
based on properties of soils as found in the field. Earlier classification 
systems were based on soil genesis rather then field properties. Utilizing 
the information developed by Stoner et. al. (1979), which includes reflectance 
of various soil samples from around the continental United States as well as 
numerous physical and chemical properties, a statistical analysis was carried 
out to determine which properties affect which wavebands. The chosen wavebands 
for this study were those which match the bands of the Thematic Mapper. The 
thematic mapper operates on six narrow bands: Band 1 (0.45-0.52?m) , Band 2 
(0.52-0.60,'tm) , Band 3 (0.63-0.69^0), Band 4 (0.76-0.90(jm) , Band 5 (1.55- 
1.75^im), Band 7 (2.08-2.35^o) , and a broad thermal band (10.4-12.5>un) . 

Since the laboratory reflectance values cover only the wavelengths 0.4 to 
2.4jnn, this study does not consider the thermal band. 

The work of Crouse, et. al. (1983) has shown that bands 5 and 7, which 
are in the middle Infra-red portion of the spectrum, are significant in the 
study of soils using satellite information. The near Infra-red band, Band 
4, also shows promise for these studies. The two visible bands, Band 2 and 
Band 3, will sometimes provide additional information. Band 1 does not 
appear to be useful for these purposes, because the laboratory data is 
sketchy in these wavelengths. For that reason Band 1 was not considered in 
the present study. 
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The statistical work in this study was carried out using the Statistical 
Analysis System package developed by SAS Institute, Cary, North Carolina. 
The procedure R-SQUARE was used to carry out the initial regression work. 
The GLM procedure was used to regress each soil property on the five TM 
bands and gave F values to indicate the significance of each TM band in the 
model. 


RESULTS 

Utilization of the published (Stoner, et. al., 1979) values limited 
this study to the following properties: Sand, Silt, Clay, Organic Matter, 
Cation Exchange Capacity, Iron Oxide content and Moisture content. By 
regressing these properties on the mean value of each of the five TM bands 
chosen for this study (Bands 2, 3* 4, 5, and 7) it was determined that 
prediction of certain properties was possible just by knowing these five 
wavebands. The results of this regression study, on the 467 soils in the 
data set , are presented in Table 1 . The R-SQUARE values can be read as 
percentages of the various properties which can be accounted for by the five 
wavebands. For instance, from Table 1 we can say that we can account for 
44.4$ of the organic matter in the soil by knowing the reflectance values of 
the five bands. 

The next step in the study was to stratify the total sample by Prent 
Material, Temperature Zone, Moisture Zone, Taxonomic Order and Climate 
(comined temperature and moisture). The results of those regressions are 
shown in Table 7 under the appropriate column headings. The values shown on 
Table 7 are mean values of the various classes within each group. Parent 
material was broken down into the following classes: Igneous Residuum, 

Sedimentary Residuum, Eolian (air transported) origin, Unconsolidated 
Terrigenous Deposits (alluvial and lacustrine), Marine Deposits, and Glacial 
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origins. Table 2 shows the various R-SQUARE values for these classes. The 
soils were then reclassified by temperature zones and analyzed. The results 
of that stratification and regression are shown in Table 3» Other classifi- 
cations based on Moisture Zone, Taxonomic Order and Climate are tabulated in 
Tables 4, 5, and 6 respectively. The summary shown in Table 7 indicates 
that prediction values increase with stratification. Looking at the values 
in Table 7 it is apparent that a stratification based on Climate (tempera- 
ture and moisture) would be beneficial. It is also intuitively apparent 
that these two factors are the only ones of which we could have prior 
knowledge for an unknown earth area. It would be easier to assign a climatic 
parameter to an unmapped area than to use Parent Material or Taxonomic 
Order, since these parameters imply a prior study of the soil. In simpler 
terms, if we knew these parameters we would no longer have an unmapped area. 

A known Taxonomic Order for a soil on the earth could, however, be used to 
verify or "fine tune" our TM data. We can assign a climatic parameter to an 
unknown earth area without any prior kne Ige of the soil type. 

The next phase of the study involved regressing the soil properties on 
the five TM wavebands in order to determine which bands are significant in 
predicting the soil properties . Table 8 shows the results of this procedure 
for the entire sample population of 467 soils. A high F value and a low 
(less than 0.0050) PR > F value indicates significance. It can be seen from 
Table 8 that bands 2 and 7 are significant in predicting six of the seven 
properties. Band 3 is significant for five properties while Band 5 and 
Bands 4 are significant to a lesser degree. The F value tests how well the 
model as a whole accounts for the dependent variable's behavior. The value 
labeled "PR > F" is the significance probability - if this value is small it 
indicates significance. 


Again, the soil sample population was stratified and the results of the 
regression analyses were tabulated. Tables 9 and 10 show the values for 
each property in relation to the five TM bands, when stratified by 
temperature (Table 9) and moisture (Table 10). These results show only the 
greatest F values, where the significance probability is les'* than .0050 for 
the various properties on the five wavebands. Table 11 shows the results of 
combining temperature zones and moisture zones for a stratification based on 
climate. The importance of certain wavebadnds to the model are obvious from 
a perusal of Table 11. Based on all stratifications, it can be seen that 
Band 2, Band 7 and Band 3 are most significant in the prediction of soil 
properties. Bands M and 5 are slightly less important overall; but they 
should be considered in all models. 

A preliminary approach to prediction of certain properties was 
attempted based on correlations which were evident when the angle of the 
tangent to the rising limb of the reflection curve was compared with the 
soil properties. This angle was designated "Angle X" and is referred to as 
such in this paper. The unstratified data indicated that correlations exist 
in Sand Content, Clay Content, Organic Matter, Cation Exchange Capacity and 

s 

Moisture Content. When the data is stratified by climate, strong 
correlations were found for C.E.C., Organic Matter, Clay Content, Sand 
Content and Iron Oxide Content for certain temperature and moisture 
combinations. 

Similar correlations were observed when looking at the declinatior 
angle of the tangent to the curve at the higher wavelengths of the data. 

This angle was designated "Angle Y" for this study. Correlations in the 
unstratified sample were strongest in Organic Matter, C.E.C., and Moisture 
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TABLE 1 

VALDES 

e>&T*M*et> ©V ^<eatecSSiN»Gr- Soil P*of>tfitTi«s 
OUtHC Piue T- N. 'BakiPs 

All. Sample's (4&r) hJ oki St£*tiphe-c> 


Pn_ope©Ty 

£•50 OAee 

5 akjp 
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51lt 

.049 

Clay 

.361 

OK&AMC 

Matter 

.444 

CA*noM 

^XCHAlOfie 

capacity 

.340 

XeoM 

.119 

OXIDE 1 

Moistubc 

CoWTiSYJfT. 

.347 
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OC6,. 
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! 

.?56 
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TA6UE 5 
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Content. The correlations were all negative-indicating that lower 
quantities of these properties produce a steeper angle of descent at the 
trailing end of the reflectance curve. When the data is stratified by 
climate, the correlations are strongest in certain climate classes for the 
three soil properties mentioned as well as a very strong relationship of 
iron oxide in the Thermic-Arid and Thermic-Semiarid sample set. Figure 2 
shows these relationships in a summarized format. 

CONCLUSIONS 

The results of this work will aid in the utilization of TM data for 
soils investigations. We can define the wavebands which will supply the 
most significant data regarding soil properties. The study has indicated 
that bands 2 and 7 are very significant in the study of soil properties and 
their influence on TM responses. Band 2 is in the visible spectrum, while 
Band 7 is one of the middle IR bands. Bands 5 and 5 are somewhat less 
significant overall but become more important when stratification by climate 
is undertaken. Band 4, the near IR band, is individually important to some 
of the properties. Bands 2 and 4 were used to determine "Angle X', the 
slope of the tangent to the ascending limb of the reflectance curve. Bands 
5 and 7 were used to obtain "Angle Y", the slope of the tangent to the 
trailing end of the curve. 

Sand Content appears to be associated with responses in the visible 
range (Bands 2 and 3) and in the middle IR bands (Bands 5 and 7). In both 
cases a high sand content produces a strong reflectance response. 
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Silt content has shown consistantly poor potential for prediction. 
Only when stratified by climate does Band 2 come into play in several 
classes to produce positive correlations. 


Clay content appears to be related negatively to responses in Bands 7, 
2, and 3. 

• Organic Matter Content, under various stratifications, figures 
dominantly in the Band 2 and Band 3 responses. These are negative 
correlations indicating that high organic content tends to depress the 

* 

responses in the visible bands. This reduced response in the visible area 
produces a flattened curve. This was used by Stoner (1981) to classify the 
curves by shape. 

Cation Exchange Capacity shows excellent results in Bands 2 and 7» 
while several stratification classes show strong responses in Bands 3 and 4. 
These responses all show negative correlations, indicating that an inverse 
relationship exists between C.E.C. and waveband response. 

Iron Oxide Content showed somewhat erratic results in the studied 
wavebands. This was not totally unexpected since Stoner pointed out iron 
absorption wavelengths in his study which lie outside the bands covered by 
the TM; however, two climatic zones (Thermic-Arid and Thermic-Semiarid) , 
showed strong negative correlations of Iron Oxide content with the Y Angle. 
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Moisture Content produced strong negative correlations with Bands 2 and 
7 with somewhat lesser results in Bands 3 and 5. Again, the negative 
correlations indicate that increasing moisture content tends to depress 
responses in these wavebands. 
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Clay content may be one key to the response prediction problem. Higher 
clay contents in soils tend to be associated with higher organic matter, 
increased Cation Exchange Capacity and increased moisture contents. The 
amount of clay present is not nearly as important as the type of clay 
minerals present. Soil scientists are quite aware of the role that climate 
plays in clay formation. All other things being equal, an increase in 
temperature and precipitation causes an increase and enhancement of the 
weathering reaction. High average temperatures tend to foster rapid 
weathering and clay formation. A minimum degree uf weathering io produced 
by climates that are: warm and dry, cold and ary, or cold and mo.ist. This 

is borne out in the climate-stratified sample population when the mean clay 
contents are examined for the various climate classes - the frigid dry, 
frigid moist and thermic dry classes all show rseu«ed clay contents. This 
indicates that it may be practical to develop a model which takes into 
account the climatic zone of an unknown earth area and then makes an assump- 
tion of clay content before checking the various wavebands for reflectance 
responses . 
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